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ABSTRACT 

A  theoretical  and  experimental  investigation  of  the  stability  of 
developing  laminar  flow  in  a  circular  pipe  subjected  to  small,  rota- 
tionally  symmetric  (torsional)  disturbances  is  presented.   Analytically, 
developing  mean  velocity  profiles  were  represented  in  polynomial  series 
form  and  subsequent  formulations  yielded  determinations  on  stability 
which  reduce  to  the  classical  conclusions  for  the  case  of  fully  developed 
flow. 

Experimentally,  torsional  disturbances  were  imposed  upon  a  developing 
laminar  pipe  flow  by  sinusoidally  rotating  a  1.5  inch  (1  diameter)  length 
of  pipe  in  the  entrance  region  of  a  carefully  constructed  wind  tunnel 
facility.   For  flows  with  maximum  Reynolds  numbers  of  8500,  maximum 
amplitude  and  frequency  of  oscillation  were  %  inch  and  30  Hz.,  respec- 
tively.  Hot  wire  anemometer  measurements  verified  the  stability  of  the 
developing  flow  fields   Moreover 9  measured  wave  speeds  and  decay  factors 
compare  favorably  to  calculated  values. 
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NOMENCLATURE 


a.   t   a.,  a_,  a_,  a.       Constants,  coefficients  of  the  power  series 

solution  of  the  disturbance  amplitude  function 

bn,  b  ,  b  ,  b  Constants,  coefficients  of  the  polynomial 

approximation  of  the  velocity  profile 

C  Dimensionless  wave  speed 

r 

C.  Amplification  factor 

d  Inside  diameter  of  the  pipe 

P  Pressure 

r  Radius  of  the  pipe 

R  Reynolds  number  (R=- — ) 

U  Mean  velocity  in  the  pipe 

U  Maximum  velocity  in  the  pipe 

max 

U  Dimensionless  velocity  profile  in  the  entry 

region,  referred  to  U 

max 

u  Axial  component  of  velocity 

w  Tangential  component  of  velocity 

z  Axial  coordinate 

°(-  Disturbance  wave  number 

l)  Kinematic  viscosity 

Stream  function 


y 
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I.   INTRODUCTION 

The  purpose  of  the  work  reported  here  was  the  study  of  the  evolution 
of  small,  rotationally  symmetric  disturbances  introduced  into  the  laminar 
transition  region  (the  hydrodynamic  entrance  length)  of  Poiseuille  pipe 
flow. 

Since  the  pioneering  work  in  the  field  of  hydrodynamic  stability 
conducted  by  0.  Reynolds  in  the  1880' s,  the  stability  of  fully  developed 
pipe  flow  has  been  the  subject  of  many  theoretical  and  experimental 
studies.   All  of  the  analyses  up  to  the  present  concerning  the  stability 
of  fully  developed  pipe  flow  show  that  such  a  flow  is  stable  to  small 
disturbances. 

A  review  of  the  theoretical  and  experimental  determinations  con- 
cerning the  transition  phenomenon  leads  one  to  the  following  considera- 
tions:  (1)  transition  impends  because  the  disturbances  introduced 
contain  some  initially  large  disturbances  so  that  small  disturbances 
are  not  significant  regardless  of  when  or  where  they  are  introduced; 
or,  (2)  initially  small  disturbances,  either  axial ly -symmetric,  rota- 
tionally symmetric  or  three-dimensional,  grow  in  size  along  a  favorable 
length  of  conduit,  namely  the  entrance  length,  via  convective  inter- 
action while  the  boundary  layer  is  rapidly  developing.   Consequently, 
these  small  disturbances,  which  would  not  have  grown  in  a  fully  developed 
region,  now  acquire  a  magnitude  which  eventually  renders  the  fully  de- 
veloped flow  unstable. 

A  theoretical  analysis  in  support  of  the  above  considerations  has 
been  presented  by  Tatsumi  [1].   His  conclusions  demonstrated  that  small 
disturbances  magnified  in  the  entrance  region  of  a  pipe  flow  field  above 
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a  minimum  Reynolds  number  of  19,400,,  do  cause  transition  after  a  passage 
of  17  tube  radii .   The  significance  of  his  study  is  not  the  determination 
of  a  value  for  the  critical  Reynolds  number,,  but  that  there  is  a  pre- 
diction of  instability,, 

Presently9  the  case  of  developing  axially-symmetric  flow  is  treated. 
Results  are  shown  to  be  reducible  to  the  determinations  of  Pekeris  [5] 
for  fully  developed  pipe  flow.   Additionally,  experimental  verification 
of  calculated  wave  speeds  and  amplification  factors  is  presented. 
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II.   THEORETICAL  ANALYSIS 

A.    CYLINDRICAL  ORR-SOMMERFELD  EQUATIONS 

To  describe  a  general  disturbance  in  mathematical  form  requires  the 
formulation  of  a  three-dimensional  (r,  0,  z)  vector  function  of  time,  t. 
The  generally  accepted  treatment  for  small  disturbances  is  to  view  them 
as  a  Fourier  series  of  a  large  number  of  elementary  disturbances,  each 
of  which  can  be  represented  in  time  and  space  by  a  finite  number  of  para- 
meters.  Indeed,  this  is  the  foundation  of  the  small  disturbance  theory  - 
that  the  disturbance  is  a  linear  combination  of  spatially  periodic 
components  which  grow  or  decay  with  time.   Such  a  component  is  character- 
ized by  the  stream  function, 


y>  =  P(r)  Exp  \lm<p  +  L$t  -L<*.2] 


(i) 


By  further  contending  that  the  imposed  disturbance  be  rotationally  sym- 
metric, then  the  0  dependence  is  eliminated,  and  the  representation 
reduces  to, 

Y  =  R*0  EKp[6oc(ct-^:)j  (2) 

Therein  OC   is  a  real  number  -  the  wave  number  of  the  disturbance  in  the 
fluid.   C  is  a  complex  number,  C  +  iC..   The  sign  of  the  imaginary  part 
of  C,  C,  determines  whether  amplification  (C„  <C  0)  or  decay  (C.  /"  0) 
occurs,  whereas  C  is  the  wave  speed  of  the  disturbance.   Experimentally, 
such  a  disturbance  is  difficult  to  generate,  since,  as  formulated,  it 
should  be  imposed  upon  the  entire  length  of  the  conduit  over  the  entire 
cross-section  of  the  conduit  simultaneously,  and  very  periodically  in  z 
as  well. 


Analytically,  the  disturbance  described  in  (2)  is  thought  of  as 
being  imposed  upon  a  steady  mean  flow  characterized  by  a  stream  func- 
tion, 


%    =■  f  ( rf*) 


(3) 


The  equations  of  motion  for  incompressible  flow  symmetric  about  the 
axis  of  a  pipe  are: 


(4b) 


-^ 


*_     *         ,      /    ^  > 


j 


^    r    -V^  +  "FTP  +   T?  (4d> 

It  is  convenient  to  introduce  the  non-dimensional  coordinates, 

and  define  R,  the  Reynolds  number  as  R  =  rj%   With  these  substituted, 
the  system  of  equations  (4)  becomes: 


Ml 


r  +       ^  +      ^  -    ^'  f  R  IX    "PJ  (5a) 


^  ,  atW  ^u^=.^/+l  7au<       (5b> 
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Defining  U  =  ?    ~r  and  V  =  - ^ — -^-  ,  it  is  noted  that  with  the  dependence 
on  0  dropped,  that  there  are  two  possible  types  of  symmetrical  distur- 
bances which  may  be  generated  and  which  are  not  mathematically  coupled. 
First,  there  exist  the  torsional  disturbances,  which  are  concentric 
with  respect  to  the  pipe  axis;  secondly,  there  exist  the  meridional 
disturbances,  which  occur  in  planes  passing  through,  and  parallel  to, 
the  axis  of  the  pipe.   Treatment  of  each  case  is  necessary  due  to  the 
aforementioned  difficulty  in  experimentally  imposing  the  mathematically 
exact  disturbance  upon  a  flow  regime.   It  is  assumed  that, 


f(r^t)     =  f(r)  Exp[d*(Ci-2)] 


(6a) 


Sl(*At)       =6COEXp[£o<(c:t-*)J  .     (6b) 

These  relations,  together  with  the  definitions  of  velocity  components 
in  terms  of  the  stream  functions,  lead  to  the  following  uncoupled  sys- 
tem of  exact  equations  of  motion: 

fcc*p+a£#  +  igg  ^D>-$  fet?V]  +  (7a) 

where      £)*  -    £.    _  £  ^_      j£_ 

As  y   is  composed  of  the  main  flow  with  the  disturbance  superimposed 
upon  it,  it  can  be  described  as. 
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r  =  %  -4-  Y 


(8) 


where    yy   refers  to  the  main  flow  and   y   is  as  defined  in  (6a). 
Substituting  (8)  into  the  system  (7),  dropping  squares  and  products 
of  (P       and  J~l^     ,  as  each  is  assumed  small ,  leads  to  the  following 
uncoupled  system  of  linearized  equations: 


ml  +.  ±  rv&  ix  ■— ""^a  iifl      =  -4-obi 


Vk        H  ^r    ~^*         *=•  .5v 


ft 


(9b) 


To  proceed  further  requires  that  the  exact  nature  of   j0     be  known. 
In  the  fully  developed  case,  j0      is  a  simple  function  of   jT  alone, 
and  is  easily  treated.   In  the  case  of  developing  flows,  such  as  the 
one  under  investigation,  there  is  no  precise,  known  functional  descrip- 
tion of0 


%    =     %(r^t) 


B.    DEVELOPING  FLOW  STABILITY  DETERMINATIONS 

Schiller  [2],  and  Langhaar  [3]  present  workable  models  of  developing 
velocity  profiles  in  parametric  form.   However,  after  long  and  fruitless 
attempts,  the  author  has  concluded  that  there  is  no  simplified  mathematical 
solution  suitable  for  representation  of  the  developing  mean  flow  at  pre- 
sent.  Accordingly,  following  Tatsumi  [1],  purely  parallel  flows  which 
have  the  same  velocity  profile  at  a  given  axial  distance  as  the  developing 
flow  shall  be  treated.   For  such  velocity  profiles,  it  was  found  that 
experimental  data  could  be  approximated  to  an  excellent  degree  of 
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accuracy  by  a  third  degree  polynomial  obtained  from  a  least-squares 
curve  fitting  procedure.   To  make  the  solution  more  general,  it  is 
assumed  that 

UCH  =   b.  +  bir  +  b^  +  l>,  **     (10) 

from  which 

^olr)    -  1ok%.  h±y\  JbaLy  +.i^-K*"     (ii) 

is  obtained.   By  obtaining  the  necessary  terms  from  (11)  and  (6)  and 
substituting  them  into  (9) s  the  following  equations  are  obtained: 

where         «  _    Cu  ^     i]    ^     ^  T  (12b) 

and  primes  denote  differentiation  with  respect  to  r  .   Treating  (12a) 
first,  it  is  clear  that  for  large  R,  the  terms  on  the  right  side  of  the 
equation  will  predominate,  since  even  when  ^C  is  assumed  small  £(R  is 
Still  large.  Yet,  if  the  entire  L? jj  term  is  ignored,  it  is  impossible 
to  satisfy  all  of  the'  boundary  conditions  of  the  problem.   Order  of 
magnitude  considerations  allow  for  the  retention  of  the  g'*  term,  and 
thus  satisfaction  of  all  boundary  conditions.   Returning  to  the  defini- 
tion of  g, 

Inserting  (13)  into  (12b)  then  yields, 

r-  *r#r-  *x- **)*'+(:$>  +-HX' (14) 


4- 

where 


and 


4£    =      ^    /"Sk^K-    ^ 
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Before  simplifications  equation  (14)  is  of  the  general  form, 

a  form  which  possesses  a  regular  singular  at  the  origin,  r  =  0.   This 
fact  guarantees  the  existence  of  at  least  one  series  solution  of  the 
form  f\  s  ■__   i^^  ^Tnu  I/-K  (15) 


Mr   K"J*«r 


£*3 

which  will  converge  for  0  <C  r  <  r  ,  where  r  is  the  distance  from  the 

o         o 

origin  to  the  next  nearest  singular  point.  As  the  disturbances  intro- 
duced are  confined  to  the  vicinity  of  the  pipe  wall,  r  sv  1,  and  again 
following  Schlichting  [6] „  the  terms 


1  f"         5.    r"       ^^ 


i/ 


and   -ys  f 


are  neglected  as  they  are  all  of  the  order  1,  whereas  the  o£R  term  is 
assumed  quite  large.   The  resultant  equation  is 

Assuming  a  solution  of  the  form  (15),  differentiating  and  inserting 


terms, 


(<r?+K   ){<*  ^c-lj^+iC-i  J[^  +  /c^  )    c^    K 


^4*-^ 


.-  o 


_<^<R 


_    *--°  J  (17) 


4,  i  *-f? 


<=o 


JT(^k)^  ^  (C-^-^-fc^-V') 


+  t*K 


O0  ■"■j   _, 


-pb»r-^)     = 


o 
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Setting  the  coefficient  of  the  highest  power  of  r  equal  to  zero  results 
in  the  indicial  equation,, 


m("i-±  )(/^-i)(/>z-3  )<*■< 


=  0 


with  indicies  n  ■  0,1,2,3  to  guarantee  that  a  :fc  0. 

Before  proceeding  further,  specific  attention  is  due  to  the  boundary 
conditions  of  the  problem.   The  no  slip  condition  at  the  wall  guarantees 
that  f'(l)  =  0,  i.e.,  — ^"Tp  =  0.   Since  the  flow  is  axially  symmetric, 
v'  must  be  symmetric,  hence  zero  on  the  axis,  or  f(0)  =  0.   Finally, 
there  can  be  no  sources  at  the  wall,  and  therefore  j   (1)  =  0,  or  f(l) 
=  0.   Here,  it  should  be  noted  that  most  previous  works  also  include 
a  condition  that  u'(0)  ■  0,  or  that  f'(0)  -  0.   However,  as  noted  by 
Leite  [4],  there  is  no  guarantee  that  only  those  disturbances  generated 
by  the  experimenter  are  being  accepted  by  the  flow,  i.e.,  there  is  no  a 
priori  reason  to  expect  u'(0)  =  0.   It  is  noteworthy  that  it  is  not  neces- 
sary to  impose  this  artificial  condition  in  this  analysis. 

Taking  the  largest  exponent  of  r  first,  it  is  found  that  the  second 
highest  ordered  coefficient, 

{n-1  JoO^t-hiynt^)^  =0  (is) 

can  only  be  identically  satisfied  of  a,  =  0.   Continuing,  for  each  co- 

1 

efficient  in  turn,  it  is  found  that  the  recurrance  relation  yields: 

1v=     c*ft  (c-boW-a\  (19) 


a^^lsLELki  frt-i)    ex*  (20) 
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(^^i)(^i-i)  <21) 


After  substituting  n  =  3,  then; 
0*  - 


fy    ^^(C-bo) 


o-x  -  -  2-  c  <  R  £, 


a4  = 


pvo 


0   ^3  ■     To 


s 


With  these  coefficients,  the  final  form  of  the  solution  becomes, 


^ 


(22) 


The  existence  of  additional  linearly  independent  solutions  is  deter- 
mined by  returning  to  the  indicial  equation.   The  value  n  =  0  will  yield 
such  a  solution,  but  as  the  boundary  condition  f(0)  =  0  would  not  be 
satisfied  without  a  =0,  then  this  solution  may  be  discarded.   The  same 
consideration  applies  for  the  value  n  =  1,  i.e.,  the  f'(l)  =  0  boundary 

condition  would  not  be  satisfied  without  a  =  0.   In  each  of  these  cases, 

o  ' 

if  a  were  to  be  equal  to  zero,  then  it  is  evident  that  the  series  could 


be  represented  as 


oo 


A\+«+i 


Therein,  the  values  n  =  0  and  n  =  1  would  not  then  satisfy  the  resultant 
indicial  equation. 

Finally,  the  case  of  n  =  2,  the  last  of  the  four  index  values  is 
considered.   Now,  n  =  2  identically  satisfies  the  recurrance  relation 
involving  a.,  i.e.,  the  second  highest  ordered  term.   Thus,  it  is 
possible  that  this  value  would  yield  a  linearly  independent  solution. 


20 


However,  as  this  exponent  differs  from  the  highest  by  a  positive  integer 
(vis.  1),  it  is  always  possible  to  generate  another  linearly  independent 
solution  of  the  form  , 


]>   CK^V  $m{r)^dKr 


aK  r 

This  solution  clearly  fails  the  f(0)  =  0  boundary  condition,  hence  it  is 
also  discarded.   Therefore,  of  the  four  linearly  independent  solutions 
possibly  generated  by  the  exponents  from  (18)  ,  only  one  remains  that 
satisfies  all  the  conditions  of  the  problem,  viz.,  the  solution  represented 
in  relation  (22). 

Applying  first  the  wall  boundary  condition  to  this  solution  yields, 

f\r)  =.  1  -  *£  cc  "t  ^§c^  -  ^$  [WrW>  #  J 
+  4y?  c-fc,  -  c^  _  t^-i^|;et-M    (23) 


Separating  this  equation  into  real  and  imaginary  parts  gives 


-     3        C^^Cc(C^bo)    ~    O 


?3       £-vo  %<+o  y.syo  "" 


/t    -  <*fflj  4    «*/*fr _    b^il  r(r 


~C    (24b) 


Solving  (24a)  for  C.  results  in 

i 


Hildebrand,  F 
Prentice-Hall,  1962 


3*£ 


(25) 


'•  »  Operational  Calculus  for  Applications ,  p.  137, 
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Following  the  same  procedure  as  above  in  applying  the  boundary  condition 
f'(l)  =  0  to  the  solution  in  relation  (22)  yields  two  additional  equa- 
tions 

Ct-b„   _  h,   _  Ja.  _  **(Ce-M_*£Qfe-Q-.g  (26a) 
»       3"    ¥0       />«  /<ro 


Solving  (26a)  for  C„  gives 

£■-  "^g  -  3g  fe/     4^>    —  J^l      r27> 

It  should  be  noted  that  C,  and  C  have  each  been  normalized  with  respect 

i      r 

to  U    and  that  ®C    has  been  normalized  with  respect  to  d.   It  is  clear 
max 

that  the  normalized  C  value,  the  wave  speed,  takes  on  values  0  ^  C  <  1 

r  r^ 

with  the  maximum  value  occurring  in  the  region  of  maximum  flow  velocity. 
Inspection  of  (25)  and  (27)  reveals  that,  within  the  developing  region, 
C,  must  always  be  positive.   Hence  the  developing  region  is  stable. 
Further,  in  the  limiting  case,  viz.,  in  the  fully  developed  flow  regime, 
the  coefficients  of  the  polynomial  representation  of  the  velocity  pro- 
file become; 

bQ  -  i,  b3  -  -l.  b2  =  b4  =  0 

Analagous  considerations  demonstrate  then  that  in  the  limiting  case  of 
fully  developed  flow,  the  solution  represented  by  (22)  also  holds. 


See  Also  reference    [5]. 
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Data  obtained  from  the  curve  fitting  procedure  is  presented  later 
to  verify  the  applicability  of  this  solution  to  the  developing  flow 
regime.   The  solution  of  (24b)  would  follow  from  Pekeris  [5]  under  the 
assumption  of  purely  parallel  flows.   Further  pursuit  of  the  exact  solu- 
tion of  (24b)  however,  depends  on  the  exact  determination  of  v  as  a 
function  of  r  and  z. 
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XII.   EXPERIMENTAL  ANALYSIS 

A.    TUNNEL  FACILITIES  AND  EXPERIMENTAL  PROCEDURE 

The  experiments  were  conducted  in  a  horizontal  Lucite  pipe  1.5 
inches  I.D.  and  45  feet  (300  diameters)  long.   A  centrigugal  fan  was 
used  to  pull  air  from  the  stagnation  conditions  in  the  plenum  chamber 
through  the  system.   The  inlet  chamber  consisted  of  a  rectangular  box 
3  feet  square  in  cross-section  and  4  feet  long  which  was  fitted  internal- 
ly with  6-24  mesh  wire  screens ,  followed  by  a  2  foot  long  settling  cham- 
ber.  Sizing  of  the  chamber  was  fixed  to  limit  the  maximum  interior 
velocity  to  1.5  fps.   The  well  rounded  inlet  horn  was  a  quarter-section 
of  a  9  inch  by  18  inch  ellipse  molded  of  fiber  glass  and  polished  smooth. 
All  joints  were  hand  fitted  to  insure  smoothness  and  uniformity  of  mating 
interior  surfaces. 

The  disturbance  generator,  Figure  1,  was  located  12  inches  (8  dia- 
meters) from  the  pipe  entrance,  and  was  followed  by  measuring  ports  at 
16,  20,  32,  40,  48,  56,  and  272  diameters  downstream  of  the  entrance. 
Uniformity  of  the  internal  diameter  and  concentricity  of  the  pipe  were 
found  to  be  within  .005  inch  for  all  sections  of  the  pipe.   The  indivi- 
dual pipe  sections  were  joined  by  flanges,  each  pair  hand  fitted  and 
joined  to  insure  uniformity.   The  pipe  sections  joined  to  the  distur- 
bance generator  were  fitted  with  rubber  0-rings  to  assist  in  damping  any 
vibrations  which  might  carry  over  from  the  driver  to  the  pipe.   The  pipe 
was  fastened  to  a  horizontal  metal  track  which  was  bolted  to  the  con- 
crete deck.   The  vibrator  was  securely  fastened  to  the  deck  after  being 
isolated  with  one  inch  of  corrugated  rubber  vibration  mounting.   The 
entire  apparatus  was  longitudinally  aligned  with  a  high  power  transit 
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and  leveled.   Maximum  deviation  from  true  horizontal  and  longitudinal 
alignment  is  felt  to  be  no  more  than  .25  inch  over  the  entire  length. 

Measurements  were  taken  with  a  Thermo -Systems  constant  temperature 
hot  wire  anemometer  utilizing  a  boundary -layer  type  probe  with  a  .095  inch 
long,  .00015  inch  diameter  tungsten  element.  The  probe  could  be  located  at 
any  of  the  above  mentioned  stations,  and  accurately  positioned  to  within 
.050  inch  of  the  pipe  wall. 

Prior  to  the  commencement  of  data-taking,  a  careful  survey  was  taken 
to  determine  the  physical  limitations  of  the  apparatus.   It  was  found 
that  the  maximum  attainable  Reynolds  number  was  9650;  however,  at  this 
setting,  the  sensitivity  of  the  flow  regime  was  such  that  any  disturbance 
within  the  laboratory  was  sufficient  to  cause  transition  of  the  flow. 
The  maximum  usable  Reynolds  number  was  found  to  be  8500.   The  survey  fur- 
ther revealed  that  above  a  frequency  of  oscillation  of  30  Hz.,  too  much 
vibration  was  being  carried  over  to  the  conduit.   Hence  data  runs  were 
limited  to  oscillations  at  30  Hz.  and  less. 

The  following  procedure  was  used  in  taking  data  during  the  experi- 
mental phase : 

1.  The  mean  velocity  corresponding  to  the  desired  Reynolds 
number  was  measured  and  monitored  at  the  station  located  in  the  fully 
developed  flow  region  (272  diameters  from  the  entrance)  by  regulating 
the  amount  of  air  pulled  by  the  fan. 

2.  The  probe  was  positioned  at  the  station  nearest  the  genera- 
tor, and  rotated  until  the  position  fixing  the  maximum  reading  was  found. 

3.  A  velocity  profile  was  first  obtained;  data  taking  runs 
then  commenced  with  oscillation  frequencies  of  30  Hz.,  with  %  and  \   inch 
amplitudes,  oscillations  then  at  20  Hz.,  with  %  and  \   inch  amplitudes, 
and  finally  at  10  Hz  at  %  inch  amplitude. 
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4.  The  probe  was  then  moved  to  the  next  location  down  stream, 
and  steps  2  and  3  repeated. 

5.  Upon  completion  of  recordings  at  all  stations  in  the  entry 
length  region  at  the  desired  Reynolds  number,  step  1  followed  by  2,  3  and 
4  were  again  repeated  for  another  Reynolds  number. 

The  following  procedures  were  used  in  recording  data: 

1.  The  output  of  the  hot  wire  was  displayed  on  an  RMS  meter 
and  on  a  dual  beam  oscilloscope,  as  was  the  signal  input  to  the  vibrator. 

2.  As  the  probe  was  moved  from  its  closest  proximity  to  the 
wall  towards  the  center  of  the  tube  in  .025  inch  increments,  the  hot 
wire  bridge  voltage  and  RMS  voltage  were  recorded  and  the  oscillogram 
either  described  or  photographed.   This  step  was  repeated  at  each  loca- 
tion for  each  data  run  conducted  in  (3)  above. 

The  data  obtained  in  the  manner  just  described  was  normalized 
through  the  use  of  appropriate  values.   Graphical  presentations  of  this 
data  are  discussed  in  the  following  sections  of  this  report. 

B.    TUNNEL  PERFORMANCE  AND  EXPERIMENTAL  RESULTS 
1.   Developing  Velocity  Profiles 

Figures  2-5  show  the  velocity  profiles  obtained  both  at  high 
and  low  Reynolds  numbers  at  two  points  in  the  developing  flow.   They  are 
compared  with  the  analytical  approximations  of  Langhaar  [2],  with  a  cor- 
rection factor  of  .8  applied  to  the  development  length.   The  correction 
factor  was  introduced  because,  similarly  to  Leite  [3],  it  was  found  that 
the  peak  velocity  obtained  at  a  given  axial  location  was  always  higher 
than  that  predicted  by  theory.   It  is  additionally  noted  that  the  boundary 
layer  itself  does  not  develop  as  rapidly  as  predicted,  although  this 
factor  is  strongly  dependent  upon  Reynolds  number  and  thus  may  be  a 
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characteristic  of  the  apparatus.   Figures  22  -  25  show  the  polynomial 
fitted  by  a  least-squares  technique,  the  results  of  which  were  used  in 
the  determination  of  the  damping  factor.   The  second  curve  in  each  of 
these  figures  is  that  of  a  fully  developed  profile  offered  for  compari- 
son.  Table  1  contains  a  listing  of  the  coefficients  of  the  polynomials 
used. 

The  noise  level  in  the  developing  flow  is  very  much  greater,  by 
approximately  an  order  of  magnitude,  than  that  previously  reported  in 
fully  developed  flow  experiments.   It  is  felt  that  due  to  the  type  of 
probe  utilized  and  because  of  it's  orientation,  some  of  the  noise  signal 
is  indeed  a  measure  of  the  radial  component  of  velocity,  v.   In  any  case, 
the  average  turbulence  level  for  all  stations  was  .8  -  1.07o.  As  might 
be  expected,  the  higher  values  occurred  at  close  proximity  to  the  wall, 
and  the  minimum  value  was  found  in  the  core  region  of  the  profile. 
The  very  nature  of  the  hydrodynamic  entry  region,  where  most  of  the 
frictional  losses  are  known  to  occur  in  pipe  flow,  suggest  that  one 
should  expect  a  rather  noisy  signal. 

2.    The  Disturbance 

At  the  highest  frequency  and  highest  amplitude,  the  turbulence 
level  jumped  to  a  peak  of  127,  near  the  wall,  with  the  peak  at  the  center 
line  of  8%  at  the  farthest  station  located  in  the  entry  length.   At  the 
highest  Reynolds  number  tested,  the  level  at  the  station  in  the  fully 
developed  region  reached  a  peak  value  of  1.07,  compared  with  .  87»  without 
any  disturbance.   It  was  noted  that  at  all  Reynolds  numbers  and  at  fre- 
quencies and  amplitudes  accepted  by  the  flow,  that  the  radial  propaga- 
tion of  the  disturbance  was  strongly  dependent  upon  the  velocity  profile, 
i.e.,  upon  the  boundary  layer  thickness.   Penetration  into  the  core  region 
was  strongly  amplitude  dependent. 
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a.  Reynolds  Number  Dependence 

One  would  expect  that,  for  a  fixed  frequency  and  amplitude, 
the  level  of  turbulence  at  any  downstream  position  along  the  pipe  would 
increase  with  increasing  Reynolds  number.   Figures  6-9  show  plots  of 
Turbulence  Level  vs.  r/a  at  near  and  far  stations  at  high  and  low 
Reynolds  numbers  and  demonstrate  this  relation  quite  well.   Figure  10 
shows  the  confining  effect  of  the  thinner  boundary  layer  found  at  higher 
Reynolds  numbers,  i.e.,   the  higher  level  disturbances  are  not  allowed 
to  diffuse  towards  the  center  into  the  core  region.   This  effect  de- 
creases dramatically  with  decreasing  Reynolds  number.   Not  only  is  diffus- 
ion allowed,  but  it  occurs  sooner,  as  the  lower  mean  velocity  would 
indicate.   Lastly  it  is  noted  that  there  were  combinations  of  frequencies 
and  amplitudes  which  produced  extremely  small  fluctuations  in  the  mean 
flow,  just  above  the  background  noise  level.   Such  phenomenon  were 
observed  at  10  Hz.  at  an  amplitude  of  %  inch  and  at  20  Hz.  at  \   inch  at 
all  Reynolds  numbers.   It  appears  that  these  disturbances  were  simply 
damped  within  the  boundary  layer  without  producing  the  characteristic 
turbulent  spike  near  the  wall. 

b.  Frequency  Dependence 

Figures  11  -  14  show  the  very  pronounced  dependence  of 
turbulence  level  upon  the  frequency  of  the  disturbance  generator.   It  is 
noted  that  as  frequency  is  lowered,  the  speed  at  which  the  wave  travels 
down  the  pipe  should  decrease.   Hence  the  station  at  which  a  particular 
level  of  turbulence  occurs  for  a  given  amplitude  of  disturbance  at  a 
given  Reynolds  number  should  be  located  further  up  stream  as  disturbance 
frequency  is  lowered.   This  is  clearly  shown  in  Figure  15.   The  idea  of 
wave  speed  varying  directly  with  the  excitation  frequency  would  seem  to 
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suggest  that,  since  the  low  frequency  disturbance  spends  more  time  within 
the  boundary  layer,  it  interacts  for  a  longer  period  with  the  strong 
viscous  forces  acting  there.   Hence  it  is  damped  more  quickly  and  com- 
pletely.  Again  the  observation  that  the  velocity  profile  itself  strongly 
effects  the  characteristics  of  the  disturbance  in  the  fluid  is  apparent. 

c.   Amplitude  Dependence 

Figures  16  -  19  show  the  increase  in  turbulence  with 
amplitude  at  low  and  high  Reynolds  numbers.   Figure  20  shows  the  strong 
dependence  of  diffusion  upon  amplitude,  i.e.,  at  higher  amplitudes,  for 
a  given  Reynolds  number  and  frequency  of  disturbance,  the  higher  the 
amplitude  of  disturbance,  the  greater  the  penetration  of  the  disturbance 
into  the  core.   Lastly,  the  fact  that  damping  occurs  at  the  frequencies 
and  amplitudes  tested  suggests  that  there  is  an  amplitude  below  which  all 
traces  of  the  disturbance  would  damp  over  a  distance  of  a  few  diameters. 
Such  a  determination  was  not  made  here,  but  its  existance  is  foreshadowed. 
3.   The  Damping  Factor 

Following  Leite  [4],  it  is  seen  that  from  the  definition  of  the 
u  component  of  velocity  and  from  (6a), 

-m;     -  Expr-j^vcc  mi        (28) 


If  a  suitable  starting  point,  z^,  is  chosen,  this  relation  becomes, 

o 


0 


sxp  r-f  "«4*l      (29) 


or, 


0o  L    Jtf(*0) 

it**) 


(30) 
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Thus, 


Log,,*  u 


(31) 


u? 


o 
(z  -  z  )/d,  the  slopes  of  these  curves  will  be, 


If  values  of  LOG  n  —7   are  plotted  vs.  the  downstream  displacement 

o 
slopes  of  these  curves 

Then,  if  either  a  value  of  wave  speed,  C  ,  or  a  value  of  wave  number,  QC  , 
can  be  determined,  then  C„  can  easily  be  computed.  As  a  check,  equations 

25  and  26  can  then  be  solved  to  yield  comparative  values.   Figure  26 

u' 
shows  a  plot  of  LOG  _  —  vs.  (z  -  z  )/d  for  data  taken  at  Reynolds 

o 

number  of  7600  at  the  maximum  frequency  and  two  amplitudes.  As  one 

would  expect,  the  slopes  increase  as  the  center  of  the  pipe  is  approached. 

Further,  at  an  r/a  of  .9,  the  larger  disturbance  is  damped  at  a  faster 

rate. 

Equations  24  were  solved  utilizing  several  experimentally  obtained 

velocity  profiles.   Results  are  presented  in  Figure  27,  which  portrays 

the  characteristic  "hairpin"  neutral  stability  curve.   Plots  are  also 

presented  for  C  values  of  .05  and  .1  which  verify  the  stability  of  the 

flows  tested. 
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IV.   CONCLUSIONS 

(1)  Flow  in  the  laminar  transition  region  of  a  pipe  of  circular 
cross-section  is  stable  to  small,  rotationally  symmetric  (torsional) 
disturbances. 

(2)  There  exists  a  strong  dependence  of  amplification  factor  and 
wave  propagation  speed  upon  the  shape  of  the  developing  velocity  profile. 

(3)  There  is  a  dependence  of  both  wave  speed,  C  ,  and  wave  number, 
•X.   upon  radial  position,  though  evaluation  of  this  relationship  will 

require  an  accurate  description  of  v(r). 

(4)  Within  the  range  of  frequencies  and  amplitudes  tested,  there 
existed  a  peak  turbulence  level  in  excess  of  10%  very  close  to  the  wall 
of  the  pipe.   However,  disturbances  of  this  magnitude  never  reached  the 
core  of  the  flow  and  transition  of  complete  turbulence  never  occurred. 

(5)  Low  frequency  disturbances  are  completely  damped  within  the 
developing  boundary  layer. 

(6)  The  experimental  results  suggest  the  existence  of  a  "critical 
layer"  wherein  the  disturbance  is  sufficiently  strengthened  to  cause 
transition  to  impend.   Without  this  sufficient  amplification,  the  dis- 
turbance  diffuses  into  the  central  flow  region,  while  continually  being 
damped. 

(7)  Small  disturbance  theory  per  se  plays  only  a  small  role  in  the 
overall  description  of  the  transition  of  pipe  flows  from  the  laminar  to 
the  turbulent  state.   Better  understanding  of  this  transition  phenomenon 
can  only  be  achieved  through  an  analysis  which  would  describe  the  spatial 
growth  of  a  finite  disturbance  rather  than  the  temporal  growth  or  decay 
of  an  infinitesimal  disturbance.   However,  even  the  latter,  if  amplified 
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by  the  flow  regime,  quickly  reaches  an  amplitude  beyond  which  lineariza- 
tion of  the  equations  of  motion  is  no  longer  possible.   Hence,  there  is 
at  present  no  mathematical  method  suitable  for  a  truely  comprehensive 
study  of  transition.   It  appears  then  that  further  experimental  studies 
are  the  methods  whereby  more  information  can  be  gathered  on  this  subject, 
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V.   COMMENTS  AND  RECOMMENDATIONS 

The  experimental  apparatus  proved  functional,  easily  maintained  and 
operated,  and  very  reliable.  It  is  recommended  that  further  pursuit  of 
this  experiment  be  preceeded  by  the  following  alterations: 

(1)  Construction  of  a  new  entrance  section  in  accordance  with 
the  design  of  Cohen  and  Ritchie  [8]  to  facilitate  flow  developing  at  very 
high  Reynolds  numbers. 

(2)  Extension  of  the  pipe  to  allow  fully  developed  flow  to  be 
achieved  at  Reynolds  number  of  order  15,000. 

(3)  Integration  of  the  air  supply  into  a  totally  closed  system 
to  remove  the  possibility  of  interference  from  air  currents  and  tempera- 
ture variations  within  the  laboratory. 

(4)  Modification  of  the  disturbance  generator  itself  for  use 
of  frequencies  above  30  Hz.  and  amplitudes  of  %"  from  the  present  mechani- 
cal design  to  an  electromagnetic  design  to  eliminate  vibration  problems. 

After  improvements  are  made  to  the  apparatus,  prior  to  the  resump- 
tion of  experiments  in  stability,  it  is  strongly  recommended  that  an 
experimental  study  of  developing  pipe  flow  itself  be  undertaken  to  deter- 
mine : 

(1)  The  developing  velocity  components  u  and  v  as  functions  of 
the  pertinent  geometric  parameters  r  and  z. 

(2)  An  analytical  or  numerical  solution  of  the  pertinent  Orr- 
Sommerfeld  equation. 

It  is  the  author's  belief  that  the  study  of  the  stability  problem 
itself  utilizing  the  small  disturbance  theory  will  ultimately  lead  to  a 
better  understanding  of  the  problem  of  beginning  of  transition  from 
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laminar  to  turbulent  flow.   Such  a  study,  however,  must  utilize  analyti- 
cal approaches  more  compatible  with  experimental  procedure  than  those 
presented  in  the  past. 
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TABLE  1 


Coefficients  of  the  polynomial  approximations  to  the  velocity  profiles  at 

z/d  =  16   and   z/d  =  56 


z/d 


Reynolds 
number 


16 

8500 

.702197 

-.752768 

2.75123 

-2.63370 

16 

5700 

.733547 

-.657861 

2.64586 

-2.7294 

56 

8500 

.823472 

-.378369 

1.89644 

-2.34190 

56 

5700 

.87101 

-.601644 

2.70048 

-3.02695 

37 


WVl 


a  Nt-tr 


38 


Oil 


Figure  2:  Developing  Velocity  Profile;  z/d  16,  R  5700 
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Figure  3:  Developing  Velocity  Profile;  z/d  56,  R  5700 
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Figure  4:  Developing  Velocity  Profile;  z/d  16,  R  8500 
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Figure  5:  Developing  Velocity  Profile;  z/d  56,  R  8500 


42 


u' 


a 


roc 


0C3 


+    30  Hz  >l/2  IW  &MP 
D     ?0  HZ  ,!/♦  /^AMF 


*****  *  ! 

Jnnn  -+■ 

BPPBPPg 


0C5 


Dli 


r/a 


Figure  6:   Radial  Distribution  of  Disturbance  Amplitude; 
z/d  16,  R  5700 
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Figure  7:  Radial  Distribution  of  Disturbance  Amplitude; 
z/d  48,  R  5700 
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Figure  S:  Radial  Distribution  of  Disturbance  Amplitude; 
z/d  16,  R  8500 
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Figure  9:  Radial  Distribution  of  Disturbance  Amplitude; 
z/d  48,  R  8500 
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Figure  10:  Radial  Distribution  of  Disturbance  Amplitude- 
z/d  56,  R  5700,  R  8500 
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Figure  11:  Radial  Distribution  of  Disturbance  Amplitude; 
z/d  16,  R  5700 
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Figure  12:   Radial  Distribution  Of  Disturbance  Amplitude; 
z/d  48,   R  5700 
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Figure  13:  Radial  Distribution  of  Disturbance  Amplitude; 
z/d  16,  R  8500 
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Figure  14:  Radial  Distribution  of  Disturbance  Amplitude; 
z/d  48,  R  8500 
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Figure  15:  Radial  Distribution  of  Disturbance  Amplitude: 
z/d  16,  R  5700 
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Figure  16:  Radial  Distribution  of  Disturbance  Amplitude- 
z/d  32,  R  5700  .  ' 
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Figure  17:  Radial  Distribution  of  Diaturbanoe  Amplitude; 
z/d  16,  R  5700 
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Figure  18:  Radial  Distribution  of  Disturbance  Amplitude; 
z/d  48,  R  5700 
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Figure  19:  Radial  Distribution  of  Disturbance  Amplitude; 
z/d  32,  R  8500 
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Figure  20:   Radial  Distribution  of  Disturbance  Amplitude; 
z/d  48,   R  8500 
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Figure  21:  Radial  Distribution  of  Disturbance  Amplitude; 
z/d  32, $  8500 
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Figure  22:  Polynomial  Approximation  to  the  Velocity  Profile; 
z/d  16,  R  5700 
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Figure  23:  Polynomial  Approximation  to  the  Velooity  Profile; 
z/d  56,  R  5700 
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Figure  24:  Polynomial  Approximation  to  the  Veloolty  Profile; 
z/d  16,  R  8500 
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Figure  25:  Polynomial  Approximation  to  the  Velocity  Profile; 
z/d  56,  R  8500 
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Figure  26:  Decay  of  Disturbance  Amplitude,  R  7600 
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Figure  27:  Neutral  Stability  Curve 
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